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Abstract — A novel system of functions (signals) on the finite 
line, called the oscillator system, is described and studied. 
Applications of this system for discrete radar and digital com- 
munication theory are explained. 

Index Terms — Weil representation, commutative subgroups, 
eigenfunctions, random behavior, deterministic construction. 



I. Introduction 

One-dimensional analog signals are complex valued func- 
tions on the real line R. In the same spirit, one-dimensional 
digital signals, also called sequences, might be considered as 
complex valued functions on the finite line Fp, i.e., the finite 
field with p elements. In both situations the parameter of the 
line is denoted by t and is referred to as time. In this work, 
we will consider digital signals only, which will be simply 
referred to as signals. The space of signals Ti. =C(Fp) is a 
Hilbert space with the Hermitian product given by 

teWp 

A central problem is to construct interesting and useful 
systems of signals. Given a system 6, there are various 
desired properties which appear in the engineering wish list. 
For example, in various situations [1], [2] one requires that 
the signals will be weakly correlated, i.e., that for every 

1(0,^)1 «1- 

This property is trivially satisfied if & is an orthonormal basis. 
Such a system cannot consist of more than dim(H) signals, 
however, for certain applications, e.g., CDMA (Code Division 
Multiple Access) [3] a larger number of signals is desired, in 
that case the orthogonality condition is relaxed. 

During the transmission process, a signal might be 
distorted in various ways. Two basic types of distortions 
are time sliift ip{t) i—s- Lr(p{t) = ip{t + r) and phase sliift 



ip{t) i-> M„,^(t) 



(y5(t), where t,w lE Fp. The first 



type appears in asynchronous communication and the second 
type is a Doppler effect due to relative velocity between the 
transmitting and receiving antennas. In conclusion, a general 
distortion is of the type (p i-^ M^jtrip, suggesting that for 
every cp (/> E 6 it is natural to require [2] the following 
stronger condition 

!(<?!), M.^Uv?) I < 1. 



Due to technical restrictions in the transmission process, 
signals are sometimes required to admit low peak-to-average 
power ratio [4], i.e., that for every ip E & with ||iy9||2 — 1 

ma.x{\p{t)\ : t e Fp} < 1. 

Finally, several schemes for digital communication require that 
the above properties will continue to hold also if we replace 
signals from & by their Fourier transform. 

In this paper we construct a novel system of (unit) signals 
So, consisting of order of p^ signals, where p is an odd 
prime, called the oscillator system. These signals constitute, 
in an appropriate formal sense, a finite analogue for the 
eigenfunctions of the harmonic oscillator in the real setting 
and, in accordance, they share many of the nice properties 
of the latter class. In particular, the system &o satisfies the 
following properties 

1) Autocorrelation (ambiguity function). For every ip E &o 
we have 

1 if (r, w) = 0, 

<^ if {t,w)^0. 



|(^,M„U^)| = 



(1) 



2) Crosscorrelation (cross-ambiguity function). For every 
(j) (fi E &0 we have 

|(0,M,„U^)| < ^, (2) 

for every t,w E Fp. 

3) Supremum. For every signal (p E Go we have 

2 

max{\p(t)\ : i e Fp} < . 

Vp 

4) Fourier invariance. For every signal ip E &o its Fourier 
transform ^ is (up to multiplication by a unitary scalar) 
also in ©o- 

Remark 1: Explicit algorithm that generate the (split) os- 
cillator system is given in the supporting text. 

The oscillator system can be extended to a much larger 
system © e, consisting of order of p^ signals if one is willing 
to compromise Properties 1 and 2 for a weaker condition. The 
extended system consists of all signals of the form M^L^(^ for 
T, w e Fp and ip E &o- It is not hard to show that ^{&e) = 

• # (So) ~ P^- As a consequence of ([T]l and (|2]i for every 
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7^ (/) e S B we have 

l(.,^)l<^. 

The characterization and construction of the oscillator sys- 
tem is representation theoretic and we devote the rest of 
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the paper to an intuitive explanation of the main underlying 
ideas. As a suggestive model example we explain first the 
construction of the well known system of chirp (Heisenberg) 
signals, deliberately taking a representation theoretic point of 
view (see [5], [2] for a more comprehensive treatment). 



II. Model example (Heisenberg system) 

Let us denote by -0 : ^ the character ^p{t) — e p . 
We consider the pair of orthonormal bases A = {5a '■ a G ¥p} 
and ~ {tp^ : a G Fp}, where — and Sa is 

the Kronecker delta function, Sa{t) — I if t ~ a and da{t) — 
if t 7^ a. 



A. Characterization of the bases A and A^ 

Let L : 7i ^ 7i be the time shift operator Lip{t) = ip{t + l). 
This operator is unitary and it induces a homomorphism of 
groups L : Fp ^ U{Ti.) given by Lr(p{t) = Lp{t + r) for any 
TGFp. 

Elements of the basis A^ are character vectors with respect 
to the action L, i.e., Lripa — i'{'^''')'^a fo'" ™y ^ Fp. In the 
same fashion, the basis A consists of character vectors with 
respect to the homomorphism M : Fp ^ [/ {Ti) given by the 
phase shift operators M^Lpit) — ■ip{wt)(p{t). 



B. The Heisenberg representation 

The homomorphisms L and M can be combined into a single 



map TT 



Fp X F, 



'p —f U (Ti.) which sends a pair (r, w) to the 
unitary operator 7r(T,a;) = ijj(^—^Tw) o Lr- The plane 
Fp X Fp is called the time-frequency plane and will be denoted 
by V. The map n is not an homomorphism since, in general, 
the operators L^- and do not commute. This deficiency 
can be corrected if we consider the group H = V x ¥p with 
multiplication given by 

(r, w, z) • (r', w', 0') = {t + t' jW + w' , z + z' + ^{tw' — t'w)). 

The map tt extends to a homomorphism tt : H ^ U{T-l) given 
by 

7r(T, W, z) = tp {—^TW + z) Mu, O Lr- 

The group H is called the Heisenberg group and the homo- 
morphism TT is called the Heisenberg representation. 



C. Maximal commutative subgroups 

The Heisenberg group is no longer commutative, however, 
it contains various commutative subgroups which can be easily 
described. To every line i C V, that pass through the origin, 
one can associate a maximal commutative subgroup Al = 
{{1,0) eV x¥p : I € L}. It will be convenient to identify the 
subgroup Al with the line L. 



D. Bases associated with lines 

Restricting the Heisenberg representation tt to a subgroup 
L yields a decomposition of the Hilbert space H into a direct 
sum of one-dimensional subspaces Ti =® Ti^' where x runs 

X 

in the set of (complex valued) characters of the group 
L. The subspace H.^ consists of vectors ip £ H such that 
tt{1)(p = xil)^- In other words, the space consists of 
common eigenvectors with respect to the commutative system 
of unitary operators {7r(Z)};g^ such that the operator tt [1) has 
eigenvalue x (/). 

Choosing a unit vector ip^ e Ti-x for every x G we 
obtain an orthonormal basis Bl ~ {/Px ■ X G -^^}- In partic- 
ular, A^ and A are recovered as the bases associated with 
the lines T ^ {(r, 0) ; r e Fp} and W = {(0, w) : w e Fp} 
respectively. For a general L the signals in Bl are certain kind 
of chirps. Concluding, we associated with every line L C V 
an orthonormal basis Bl, and overall we constructed a system 
of signals consisting of a union of orthonormal bases 

eH={'peBL:LciV}. 

For obvious reasons, the system 6h will be called the 
Heisenberg system. 



E. Properties of the Heisenberg system 

It will be convenient to introduce the following general 
notion. Given two signals (f),(p £ Ti., their matrix coefficient is 
the function m^^^ : H ^ C given by m^,^{h) = {(j), 7T{h)ip). 
In coordinates, if we write h — {t,w,z) then m^p^^ih) = 
-0 (— irw + z) (0, Mu, o L^-iys). When (p — cp the function 



is called the ambiguity function of the vector ip and 
is denoted by A^ = m^_^. 

The system &h consists of p + 1 orthonormal basefl alto- 
gether p (p + 1) signals and it satisfies the following properties 
[5], [2] 

1) Autocorrelation. For every signal ip E Bl the function 
\A^\ is the characteristic function of the line L, i.e.. 



iv)\ = 



0, v^L, 

1, V e L. 



2) Crosscorrelation. For every (f) E Bl and p E Bm where 
L 7^ M we have 

Vp 

for every v G V. If L = M then m^^^ is the 
characteristic function of some translation of the line 
L. 

3) Supremum. A signal ip e &h is a unimodular function, 
i.e., \'p{t)\ — for every t e Fp, in particular we have 

max{|i^(t)| : i e F„} = — < 1. 

Remark 2: Note the main differences between the Heisen- 
berg and the oscillator systems. The oscillator system consists 

'Note that p + 1 is the number of hues in V . 
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of order of signals, while the Heisenberg system consists of 
order of signals. Signals in the oscillator system admits an 
ambiguity function concentrated at G V (thumbtack pattern) 
while signals in the Heisenberg system admits ambiguity 
function concentrated on a line. 

III. The oscillator system 

Reflecting back on the Heisenberg system we see that each 
vector (y9 G is characterized in terms of action of the 
additive group Ga — Fp. Roughly, in comparison, each vector 
in the oscillator system is characterized in terms of action of 
the multiplicative group G,„ = . Our next goal is to explain 
the last assertion. We begin by giving a model example. 

Given a multiplicative characte]0 x ■ Gm ^ C^, we define 
a vector x e by 

We consider the system Bstd — {x'- X'^ G^, x 7^ l}, where 
is the dual group of characters. 



C. Systems associated with maximal (split) tori 

The group Sp consists of various types of commutative 
subgroups. We will be interested in maximal diagonalizable 
commutative subgroups. A subgroup of this type is called 
maximal split torus. The standard example is the subgroup 
consisting of all diagonal matrices 

which is called the standard torus. The restriction of the Weil 
representation to a split torus T C Sp yields a decomposition 
of the Hilbert space H into a direct sum of character spaces 
^ = ® where x runs in the set of characters T^. Choos- 
ing a unit vector ip^ e H.^ for every x we obtain a collection 
of orthonormal vectors Bt — {(fy, : x & ^^^i X 7^ o'}- Over- 
all, we constructed a system 

6^={(p eBr-TcSp split}, 

which will be referred to as the split oscillator system. We note 
that our initial system Bstd is recovered as Bstd ~ Ba- 



A. Characterizing the system Bstd 

For each element a G Gm let : H ^ H be the unitary 
operator acting by scaling Paf{t) = ip{at). This collection of 
operators form a homomorphism p : Gm — * U{T-l). 

Elements of Bstd are character vectors with respect to p, i.e., 
the vector x satisfies (x) = x{o)x for every a G Gm- In 
more conceptual terms, the action p yields a decomposition 
of the Hilbert space Ti, into character spaces H = ® Ti^' 
where x runs in the group G^. The system Bstd consists of 
a representative unit vector for each space Ti,^, x 7^ 1- 

B. The Weil representation 

We would like to generalize the system Bstd in a similar 
fashion like we generalized the bases A and in the 
Heisenberg setting. In order to do this we need to introduce 
several auxiliary operators. 

Let p^ : H H, a G , be the operators acting by 
Pafit) — cr{a)(p{a~^t) (scaling), where a is the unique 
quadratic character of , let py : 7i — > 7i be the operator 
acting by prp(p{t) = %j){t^)(.p{t) (quadratic modulation), and 
finally let pg ■.H—^'H be the operator of Fourier transform 

where v is a normalization constant [6]. The operators p^, pj, 
and pg are unitary. Let us consider the subgroup of unitary 
operators generated by Pa-.Ps ^nd prp. This group turns 
out to be isomorphic to the finite group Sp — SL2{¥p), 
therefore we obtained a homomorphism p : Sp U{TL). 
The representation p is called the Weil representation [7] and 
it will play a prominent role in this paper 

multiplicative character is a function x '■ Cm — > which satisfies 

X(^y) = X(^)x(y) for every x,y & Gm- 



D. Systems associated with maximal ( non-split) tori 

From the point of view of this paper, the most interesting 
maximal commutative subgroups in Sp are those which are 
diagonalizable over an extension field rather than over the 
base field ¥p. A subgroup of this type is called maximal non- 
split torus. It might be suggestive to first explain the analogue 
notion in the more familiar setting of the field K. Here, the 
standard example of a maximal non-split torus is the circle 
group 50(2) C 5L2(K). Indeed, it is a maximal commutative 
subgroup which becomes diagonalizable when considered over 
the extension field C of complex numbers. 

The above analogy suggests a way to construct examples of 
maximal non-split tori in the finite field setting as well. Let us 
assume for simplicity that —1 does not admit a square root in 
Fp. The group Sp acts naturally on the plane V" = Fp x Fp. 
Consider the symmetric bilinear form _B on given by 

B{{t,w),{t',w')) = tt' + ww'. 

An example of maximal non-split torus is the subgroup 
Tns C Sp consisting of all elements g G Sp preserving the 
form B, i.e., g G Tns if and only if B{gu,gv) = B{u,v) 
for every u,v G V. In the same fashion like in the split 
case, restricting the Weil representation to a non-split torus 
T yields a decomposition into character spaces H 
Choosing a unit vector ip^ G Ti-x for every x € T"^ we obtain 
an orthonormal basis Bt- Overall, we constructed a system of 
signals 

<So'={'P ^ 13t --T c Sp non-split} . 

The system &q will be referred to as the non-split oscillator 
system. The construction of the system 60 = 6q U ©q* 
together with the formulation of some of its properties are the 
main contribution of this paper. 
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E. Behavior under Fourier transform 

The oscillator system is closed under the operation of 
Fourier transform, i.e., for every ip e &o we have ^ € So- 
The Fourier transform on the space C (Fp) appears as a 
specific operator p (w) in the Weil representation, where 

^=(-1 J)^^^^- 

Given a signal G Bt C ©o^ its Fourier transform ^ = 
p (w) ip is, up to a unitary scalar, a signal in Bt> where T' = 
wTw~^ . In fact, &o is closed under all the operators in the 
Weil representation! Indeed, given an element g G Sp and a 
signal G Bt have, up to a unitary scalar, that p (<?) 
e Bt', where T' = gTg'^. 

In addition, the Weyl element w is an element in some 
maximal torus Tw (the split type of Tw depends on the 
characteristic p of the field) and as a result signals tp G Bt„ 
are, in particular, eigenvectors of the Fourier transform. As 
a consequences a signal (p G Bt„ and its Fourier transform 
(p differ by a unitary constant, therefore are practically the 
"same" for all essential matters. 

These properties might be relevant for applications to 
OFDM (Orthogonal Frequency Division Multiplexing) [8] 
where one requires good properties both from the signal and 
its Fourier transform. 

F. Relation to the harmonic oscillator 

Here we give the explanation why functions in the non- 
split oscillator system &^ constitute a finite analogue of the 
eigenfunctions of the harmonic oscillator in the real setting. 
The Weil representation establishes the dictionary between 
these two, seemingly, unrelated objects. The argument works 
as follows. 

The one-dimensional harmonic oscillator is given by the 

differential operator D = d"^ — f^. The operator D can be 
exponentiated to give a unitary representation of the circle 
group p : SO {2,R) — > U {L^{R)) where p{e) = e'^^. 
Eigenfunctions of D are naturally identified with character 
vectors with respect to p. The crucial point is that p is 
the restriction of the Weil representation of SL2 (M) to the 
maximal non-split torus 5*0 (2,R) C SL2 (M). 

Summarizing, the eigenfunctions of the harmonic oscillator 
and functions in 6q^ are governed by the same mechanism, 
namely both are character vectors with respect to the restric- 
tion of the Weil representation to a maximal non-split torus in 
SL^. The only difference appears to be the field of definition, 
which for the harmonic oscillator is the reals and for the 
oscillator functions is the finite field. 

IV. Applications 

Two applications of the oscillator system will be described. 
The first application is to the theory of discrete radar. The 
second application is to CDMA systems. We will give a brief 
explanation of these problems, while emphasizing the relation 
to the Heisenberg representation. 



A. Discrete Radar 

The theory of discrete radar is closely related [2] to the 
finite Heisenberg group H. A radar sends a signal (p{t) and 
obtains an echo e{t). The goal [9] is to reconstruct, in maximal 
accuracy, the target range and velocity. The signal p{t) and 
the echo e{t) are, principally, related by the transformation 

where the time shift r encodes the distance of the target from 
the radar and the phase shift encodes the velocity of the target. 
Equivalently saying, the transmitted signal (p and the received 
echo e are related by an action of an element ho G H, i.e., 
e = Tr{ho)(p. The problem of discrete radar can be described 
as follows. Given a signal (p and an echo e = 'iT{ho)<p extract 
the value of ho- 

It is easy to show that |m<^_e (^)| = l^ip {h • ho)\ and it 
obtains its maximum at Hq^. This suggests that a desired 
signal (p for discrete radar should adnoit an ambiguity function 
which is highly concentrated around G H, which is a 
property satisfied by signals in the oscillator system (Property 
2). 

Remark 3: It should be noted that the system &o is "large" 
consisting of aproximately p^ signals. This property becomes 
important in a jamming scenario. 

B. Code Division Multiple Access (CDMA) 
We are considering the following setting. 

• There exists a collection of users i E I, each holding a 
bit of information 6j G C (usually bi is taken to be an 
iV'th root of unity). 

• Each user transmits his bit of information, say, to a central 
antenna. In order to do that, he multipUes his bit bi by a 
private signal ip^ gH and forms a message Ui = biPi. 

• The transmission is carried through a single channel (for 
example in the case of cellular communication the chan- 
nel is the atmosphere), therefore the message received by 
the anterma is the sum 

i 

The main problem [3J is to extract the individual bits bi from 
the message u. The bit bi can be estimated by calculating the 
inner product 

{'Pi,u) = J2 {fi, Uj) = E bj {<Pi, 'Pj) = bi+J2 bj {'Pi, p>j) . 

The last expression above should be considered as a sum of 
the information bit bi and an additional noise caused by the 
interference of the other messages. This is the standard sce- 
nario also called the Synchronous scenario. In practice, more 
complicated scenarios appear, e.g., asynchronous scenario - 
in which each message Ui is allowed to acquire an arbitrary 
time shift Ui{t) 1-^ Ui{t + Ti), phase shift scenario - in which 
each message Ui is allowed to acquire an arbitrary phase shift 
Ui{t) 1-^ e~^^^Ui{t) and probably also a combination of the 
two where each message Ui is allowed to acquire an arbitrary 
distortion of the form Ui{t) 1— > e p ^' Ui{t + Ti). 
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The previous discussion suggests that what we are seeking 
for is a large system 6 of signals which will enable a reliable 
extraction of each bit bi for as many users transmitting through 
the channel simultaneously. 

Definition 4 (Stability conditions): Two unit signals (f) ^ 
ip are called stably cross-correlated if \mip^^{v)\ <C 1 for 
every v £ V. A unit signal (p is called stably autocorrelated 
if \A^ {v)\ <C 1, for every v 0. A system 6 of signals 
is called a stable system if every signal p E 6 is stably 
autocorrelated and any two different signals <j),tp E & are 
stably cross-correlated. 

Formally what we require for CDMA is a stable system &. 
Let us explain why this corresponds to a reasonable solution 
to our problem. At a certain time t the antenna receives a 
message 

leJ 

which is transmitted from a subset of users J C /. Each 
message Uj, i £ J, is of the form Ui = 6^6 p (p^{t + Ti) — 
biT^{hi)(f^, where hi G H. In order to extract the bit bi we 
compute the matrix coefficient 

m^i.u = biRh.A^^ + #( J - {i})o{l), 

where Rh^ is the operator of right translation Rh^Aip. (h) = 
A^.{hhi). 

If the cardinality of the set J is not too big then by 
evaluating m^.^u at ft, = h^^ we can reconstruct the bit 6j. It 
follows from ([T]i and dU that the oscillator system &o can 
support order of users, enabling reliable reconstruction 
when order of ^ users are transmitting simultaneously. 

Remark about field extensions. All the results in this paper 
were stated for the basic finite field Fp for the reason of 
making the terminology more accessible. However, they are 
valid for any field extension of the form ¥q with q = p^. 
Complete proofs appear in [6]. 
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